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ABSTRACT 

We present an analysis of NGC2204-S892 - a new detached eclipsing binary composed of 
two late K dwarfs. Based on three photometric campaigns launched in 2008 we obtained 5 light 
curves (3 in V, lin B and 1 in /), and derived an orbital period of 0.451780±0. 000001 d. We 
also obtained 20 VLT/UVES spectra, enabling accurate radial velocity measurements. The 
derived masses and radii of the components (mi =0.733±0.005 Mq and _Ri =0.72it0.01 Rq; 
m2 = 0.662±0.005 Mq and R2 = 0.68it0.02 Rq) are consistent with the empirical mass-radius 
relationship established recently for lower main sequence stars in binary systems; in particular 
we find that both stars are oversized compared to theoretical models. NGC2204-S892 is very 
active: both components show variable emission in Ha and H/3 and are heavily spotted, 
causing the light curve to show appreciable changes on a timescale of weeks. Our results add 
to the increasing evidence that the observed inflation of the radii of K and M stars is related 
to high levels of magnetic activity. 

Key words: binaries: eclipsing - stars: individual: NGC2204-S892 - stars: 
low-mass - stars: K-type 

1 Introduction 

Significant discrepancies are known to exist on the lower main sequence (LMS) 
between theoretical predictions and the actual parameters of stars. For stars in 
the 0.2-0.8 Mq mass range observed radii are larger by up to 10 - 15 percent, 
and effective temperatures lower by several percent than predicted by theory. 
According to Morales, Ribas, & Jordi (2008), "... current stellar structure and 
evolution models are not adequate to describe the physical (radii) and radia- 
tive (effective temperature, color indices) properties of... low-mass stars. The 
problem is particularly severe for the analysis of open clusters or star forming 
regions". Indeed, LMS stars (LMSS) appear 50 - 90 percent older or younger 
in mass-radius diagrams than they really are, depending on whether post- or 
pre-MS models are used (Stassun et al. 2009). No definitive explanation for the 
inflated radii of LMSS has yet been given, but mounting evidence suggests that 
they result from intense magnetic activity associated with high spot-coverage 

*Based on photometric observations performed at CTIO, LCO and SAAO, and spectro- 
scopic data collected with the Very Large Telescope at ESO Paranal Observatory under pro- 
gramme 382.D-0439{A). 
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(Chabricr et al. 2007; Ribas ct al. 2008; Morales ct al. 2008). A positive cor- 
relation between radius and activity level was found for a few members of close 
binaries by Lopez-Morales (2007). This important discovery clearly deserves 
further study, extending it onto a larger region of parameter space. 

Active low-mass stars are fast rotators with rotation periods shorter than 
~3 days, and their activity is strongly correlated with rotational velocity [e.g., 
Mohanty & Basri 2003). It is evident, then, that the best objects with which 
to study the dependence of radius on activity are members of close eclipsing 
binary systems. First, their masses and radii can be accurately determined from 
photometric and spectroscopic data. Second, they are rotationally synchronized, 
and due to the spread in orbital periods they should exhibit a wide range of 
activity. Unfortunately, LMS binaries are very rare: according to Shaw & Lopez- 
Morales (2007), there is only one such system per about million stars listed in 
recent sky surveys; whereas Ribas et al. (2008) mention only 7 binaries for 
which masses and radii have been derived with uncertainties below 3 percent. 
The situation may improve owing to future large- area synoptic surveys like SDSS 
II, Pan-STARRS and LSST (Blake et al. 2008), but at the moment every new 
LMS binary with components in the mass range 0.2 M© <m< 0.8 Mq is highly 
valuable as a potential source of data enabling detailed studies of the effects of 
magnetic activity on stellar structure and evolution. 

It should be mentioned here that optical interferometry has recently enabled 
accurate direct measurements of the radii of nearby LMSS. However, the masses 
are not directly measurable in this case, and can only be estimated through 
comparison with empirical mass-luminosity relations or other indirect methods. 
Moreover, interfcrometric measiircmcnts are limited to single stars whose rota- 
tional velocities rarely exceed 3 km s^^, implying a low activity or lack thereof. 
Not surprisingly, the radii of single M dwarfs measured by interferometry do 
not show obvious excesses, and within observational errors they are compatible 
with models (Demory et al. 2009). Those authors detected possibly meaningful 
deviations only for early K stars, and found them to disappear when the mixing 
length parameter Imi^/Hp of 1.5-1.9 is used instead of 1.0. 

In the present paper we analyze NGC2204-S892 - a newly discovered de- 
tached double-lined eclipsing binary composed of two late K dwarfs (hereafter: 
S892). The binary is a foreground object in the field of the open cluster 
NGC 2204, originally described by Rozyczka et al. (2007) as star #892. They 
reported orbital period of slightly less than 0.5 d with brightness varying be- 
tween V K n A mag at maximum and V ^ 18.4 mag at primary minimum, and 
secondary minimum ~0.55 mag deep. Based on the dereddened B — V color and 
the difference between the depths of the minima, they estimated -^0.6 Mq for 
the mass of the primary. 

In 2008 we collected 5 light curves of S892 (3 in V , I in B and 1 in 
7) and derived an orbital period of 0.451780±0.000001 d. We also obtained 
20 UVES spectra enabling accurate measurements of the radial velocity ampli- 
tudes (145±2 km s^^ and 165±2 km s~^ for primary and secondary, respec- 
tively). The details of our observing runs are provided in Section 2. Section 3 
presents the analysis of the collected data and the resulting physical parameters 
of the binary. The properties of the binary are discussed in Section 4 in relation 
to the discrepancies between empirical and theoretical mass-radius relations for 
LMSS, and a brief summary of our work is presented in Section 5. 
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2 Observations 

S892 is located in Canis Major at asooo = 06^15™55?42 and (Jsooo = -18°44'51"7. 
We observed it for a week in January 2008 with the 2.5-ni Irenee du Pont tele- 
scope at Las Canipanas Observatory (LCO), for another week in November 
2008 with the 1.0-m Yale telescope at Cerro Tololo Interamerican Observatory 
(CTIO), and for 2 weeks in December 2008 with the 1.0-m Elisabeth telescope 
at the South African Astronomical Observatory (SAAO). The observations were 
made in Johnson B and V, and Cousins / bands, yielding a total of 1282 frames. 
A journal of the observations is given in Table [1] together with numbers of mea- 
surements obtained in each run and filter. The spectra of the system were 
taken between early October 2008 and mid- January 2009 with the UVES spec- 
trograph at VLT UT2 (Kueyen) at the ESO Paranal Observatory under pro- 
gramme 382.D-0439(A). 



Table 1: Journal of photometric observations 



Site 


Observing period 


Filters 


No. of frames 


LCO (2.5m) 


11.01-17.01.2008 


B/V 


85/237 


CTIO (1.0m) 


05.11-13.11.2008 


V 


464 


SAAO (1.0m) 


03.12-16.12.2008 


V/I 


331/193 



2.1 Photometry 

2.1.1 LCO 

Images were taken in the B and V bands with the TEK5 camera on the 2.5-m 
Irenee du Pont telescope at a scale of 0.259 arcsec/pixel. The FoV was 8.8x8.8 
arcmin^, and exposure times ranged from 50 to 120 s in 1/ and from 120 to 180 s 
in B. In a total of ^17.5 hours of monitoring time 85 frames in B and 237 in 
V were collected. All B frames and 235 of the V frames were useful. The PSF 
had a FWHM ranging from 0.7 to 2.3 arcsec in B and from 0.6 to 2.2 arcsec in 
V ^ with median values of 1.1 and 1.0 arcsec, respectively. 

2.1.2 CTIO 

Images were taken with the Y4K camera on the Yale 1.0-m telescope in V band 
only, at a scale of 0.289 arcsec/pixel. The FoV was 19.7x19.7 arcmin^, and 
every frame was exposed for 100 s. In a total of ^^41 hours of monitoring time 
464 frames were collected, 459 of which were useful. Only the FoV quadrant 
containing the binary was processed through the reduction pipeline. In this 
quadrant the PSF had a FWHM ranging from 1.1 to 2.2 arcsec, with a median 
of 1.6 arcsec. 

2.1.3 SAAO 

Images were taken in the V and / bands with the STE4 camera on the 1.0- 
m Elisabeth telescope at a scale of 0.31 arcsec/pixel. The FoV was 5.3x5.3 
arcmin^, and exposure times ranged from 100 to 200 s in and from 70 to 
150 s in I. In a total of ^-^37 hours of monitoring time 331 frames in V and 
193 in / were collected, out of which 322 and 181 were useful, respectively. The 
seeing was generally poor, and the resulting PSF had a FWHM ranging from 
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1.3 to 3.6 arcsec in V and from 1.2 to 3.3 arcsec in /, with median values of 2.3 
and 1.9 arcsec, respectively. 

2.1.4 Data reduction 

The preliminary processing of the raw data was performed under IRApH. All 
frames were de-biased and flat-fielded with median-averaged sky flats. The 
photometry was performed with the DAOPHOT/ALLSTAR package (Stetson 
1987). The reduction procedure started from the identification of stars with the 
subroutine FIND, followed by aperture photometry with the subroutine PHOT. 
A PSF varying quadratically with {x,y) coordinates was then constructed for 
each frame based on 30-40 isolated stars, and used for profile photometry with 
the ALLSTAR subroutine. The images were inspected visually, and for each 
camera and each filter one of the best frames was chosen as a template. Instru- 
mental magnitudes of template stars were transformed to the standard system 
using the data of Kassis et al. (1997) available in the WEBDA database. For 
each frame taken with the same camera and filter the template stars were identi- 
fied by means of transformations of (x,y) coordinates, and ALLSTAR photome- 
try for those frames was transformed to the standard system with the template 
serving as an intermediary. The quality of the transformation was verified by 
plotting the light curve of a comparison star with similar colors which was 
located close to S892 in the frames: it did not show any secular trend or discon- 
tinuity. The mean l^-band magnitude of that comparison star calculated from 
all 2008 observations was Vcomp = 17.627 mag. The rms deviation from Fcomp 
was 0.006, 0.017 and 0.016 mag for LCO, CTIO and SAAO data, respectively, 
and in all cases it was practically equal to the mean error of the photometry 
calculated by DAOPHOT. The corresponding mean magnitudes obtained from 
the three sets of data differed from T4omp by less than 10"'^ mag. 

Based on all 2008 data in the F-band we derived the following ephemeris of 
the system: 

<o = HJD 2454480.6326 ± 0.0001 (1) 
P = 0.451780 ± 0.000001 d 

The resulting phased y-band light curves are shown in Fig. [TJ Their variability 
can be immediately seen by comparing the LCO curve with any of the remaining 
two. Even the CTIO and SAAO curves differ noticeably, although they are 
separated by just one month in time. 

2.2 Spectroscopy 

UVES spectra of S892 were obtained in service mode on 10 nights spaced irreg- 
ularly between 03.10.2008 and 18.01.2009. We used the red arm of the spectro- 
graph in a standard setting centered at 5800Awith the detectors binned 2x2. 
S892 is rather faint, but its short orbital period sets stringent limits on the ex- 
posure time (long exposures would blur all spectral features). To keep both S/N 
and effective resolution as high as possible we chose a toxp of 1200 s (i.e., 3% 
of the orbital period). In practice two spectra exposed for 600 s each were 
taken one after another, and combined during data reduction. Such a proce- 
dure yielded an average S/iVsa 18 in the reduced spectra, and also facilitated 
the removal of cosmic rays. 

tiRAF is distributed by the National Optical Astronomy Observatories, which are operated 
by the AURA, Inc., under cooperative agreement with the NSF. 
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Figure 1: Phased V-band light curves obtained at (from top to bottom): LCO, 
CTIO and SAAO. The errors, with mean values of 0.007, 0.017 and 0.014 mag, 
respectively, are calculated by DAOPHOT. The zero-point on the vertical axis 
corresponds to 17.5 mag, and the curves are displaced from one another by 
0.4 mag. 

On every observing night four 600 s exposures of S892 were taken in the 
following sequence: exposure, ThAr calibration spectrum, exposure, flat-field 
image, exposure, ThAr calibration spectrum, exposure, flat-field image. Each 
spectrum was delivered in two segments recorded by independent CCD's: the 
"blue" one extending from 4780 to 5755A and the "red" one extending from 
5835 to 68OOA. We obtained 20 blue and 20 red spectra altogether, at orbital 
phases between 0.18 and 0.97, with most of them grouped around quadratures. 
The reduction of the data was performed under IRAF. The raw frames were 
debiased and flatfielded, and the spectra were extracted and callibrated with 
the help of NOAO, IMRED and ECHELLE packages. 

Radial velocity measurements were based on the broadening function (BE) 
formalism discussed by Rucinski (2002). The main cause of line broadening in 
S892 is fast rotation of the components, contributions from other factors are 
practically negligible. Eor a rigidly rotating spherical star we have 

BFrotiv) = A[{l-P)Vl-a^ + jl3{l-a^)] (2) 



with a 



V - Vrad 

Vrot sinj 



where BFrot is the rotational BE, Vrad and Vrot are, respectively, radial and 
rotational velocity, /3 is the coefHcient of linear limb darkening, and A is an 
amplitude. The observed BEs were extracted separately from blue and red seg- 
ments of all spectra, using the spectrum of K5V star HD 10361 from the UVES 
library as a template. Model BEs used for our analysis were convolutions of two 
rotational BEs given by equation ([2]) and a Gaussian with a standard deviation 
of 40 km s^^. Eits to the observed BEs were performed by means of nonlinear 
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Table 2: Radial velocities. 

HJD-2454000 Phase [km s"-^] ua [km s"-^] 



742.8072 


0.3146 


-131.6±2.0 


149.4 ± 


0.9 


742.8239 


0.3518 


-117.1±4.2 


132. 7± 


6.7 


744.7988 


0.7231 


145.6±0.4 - 


-155.5± 


5.5 


744.8155 


0.7600 


149.0±2.6 - 


-157.7± 


4.4 


746.8112 


0.1774 


-128.4±5.7 


145.7 ± 


6.3 


746.8279 


0.2144 


-141.5±2.1 


159.5 ± 


0.4 


773.7041 


0.7041 


140.4±1.4 - 


-151.9±10.9 


773.7167 


0.7319 


148.5 ±3.6 - 


-156.1± 


9.9 


775.7511 


0.2351 


-141.8±1.4 


158.6± 


5.2 


775.7679 


0.2722 


-142.6±3.5 


159.2 ± 


5.0 


792.7760 


0.9191 


73.6±0.2 


-73.0 ± 


1.8 


792.7928 


0.9562 


48.3±2.8 


-54.8 ± 


5.2 


793.8002 


0.1860 


-132.9±4.4 


151.6± 


4.5 


793.8169 


0.2230 


-142.2±3.9 


160.9± 


5.6 


796.7807 


0.7832 


146.6±6.8 - 


-161.8± 


5.3 


796.7974 


0.8202 


138.9±0.5 - 


-145. 5± 


4.9 


798.6566 


0.9354 


60.4±1.6 


-66.3 ± 


1.7 


798.6737 


0.9734 


40.6±0.9 


-30.7± 


7.8 


849.5991 


0.6950 


140. 6± 1.4 - 


-151. 2± 


8.9 


849.6159 


0.7323 


148.9 ±0.1 - 


-160.0±12.6 



Barycentric radial velocities of S892 measured from the UVES data and phased according to 
the ephemeris given by equation JTJ. 

least-squares procedures from the GNU Scientific Library. The resulting radial 
velocities listed in Tableware averages from fbiuc and Wrcd, i.e., fits to blue 
and red segment of the spectrum (the errors are estimated from the formula 
f^v = l^^biuc — i'rcd|/'\/2). As a by-product we obtained projected rotational ve- 
locities of both components of the S892 system: ui sini = 79.2± 5.9 km s^^ for 
the primary and i;2sin« = 76.2 ±8.8 km s^^ for the secondary. We note that the 
accuracy of the velocity measurements is degraded by noticeable distortions of 
the observed BFs due to stellar spots. Not surprisingly, the observational errors 
are significantly and systematically larger for the secondary component whose 
lines are weaker and more difficult to measure. 

3 Modeling 

Wc modeled the S892 system using the PHOEBE interface (Prsa & Zwitter 
2005) to the Wilson-Devinney code (Wilson & Devinney 1971). Because the 
narrow eclipses in Fig. [l]indicate that the components of the binary are far from 
contact, a detached configuration was adopted for all calculations. However, 
since the orbital period of S892 is rather short, and its light curves suggest that 
both components are nonspherical, we also included proximity effects for the 
primary and secondary. The bolometric albedo was kept fixed at 0.5, which 
is the standard value for stars with convective envelopes, and we enabled the 
double reflection option. For the limb darkening a logarithmic law was used 
as implemented in PHOEBE 0.31a. In all calculations the gravity brightening 
coefficient /3i was set to 0.32 - the classical value obtained by Lucy (1967). To 
check the sensitivity of the results to that parameter, a few converged fits were 
re-run with /3i = 0.2 - a value obtained by Claret (2000). Changes, if any, were 
of the same order as the accuracy of the fit. 
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As there was no indication for nonzero eccentricity in any of our five light 
curves, we set e = 0. The period was adopted from the photometric data. Both 
components were assumed to be entirely synchronized, ie., their synchronicity 
parameters were set to 1.0. Having no clues about the chemical composition of 
S892 we assumed that it has solar metallicity. Rozyczka et al. (2007) estimated 
that the system is located not more than ~1 kpc away from the Sun and ~250 pc 
below the Galactic plane, so that it most probably belongs to the thin disk 
population, and its true metallicity should not differ much from the assumed 
value. We re-ran a few converged fits with [Fe/H] — -0.5 and [Fe/H] = 0.5. No 
significant changes were found compared to the case with [Fe/H] = 0. 

The estimates of Rozyczka et al. (2007) were used as the initial guess for the 
masses of the components. Based on masses and period, we obtained an estimate 
of the semimajor axis a of the system. An initial guess for the inclination i was 
provided by a preliminary fit to the SAAO F-band curve, which was the most 
symmetric one among the three. A preliminary fit to the velocity curve in turn 
yielded updated (i.e., . compatible with i) approximations for masses and o. 
We checked that keeping Ti at a value fixed between 4100 K and 4300 K does 
not significantly influence the outcome of iterations for any parameter of the 
system except T2, and based on our approximate value of mi =0.75 Mq we 
adopted Ti =4200 K - in agreement with T2 =4220 ± 150 K obtained by Torres 
& Ribas (2002) for the secondary component of V818 Tau with TO2 = 0.76 or 
T2 =4200 ±200 K obtained by Bayless & Orosz (2006) for both components of 
2MASS J05162881-h2607387 with mi =0.79 and m2 =0.77 Mq. 

Because of the evident variability of the light curve (see Fig. [T]) we could not 
use the photometric data en bloc. Instead, we adopted the following procedure: 

1. Choose one of LCO, CTIO or SAAO photometric data sets. 

2. Keeping a and mass ratio q=m2/mi fixed, make a fit to the light curve by 
iterating for T2, i, primary's luminosity Li and surface potentials of both 
componets. Add spots to account for the asymmetries of light curves. The 
result will serve as the base fit for the next steps of the procedure. 

3. Fix all parameters of the base fit except a and q. Make a fit to the velocity 
curve by iterating for a and q. 

4. Go to (2). Repeat steps (2) - (4) until the base fit converges upon com- 
pleting (2). We consider the fit converged when stellar radii Ri and i?2 
change by less than 0.01 Rq during one iteration step. 

5. Fix all parameters of the converged base fit. Substitute the first photomet- 
ric data set with one of the remaining two. Adjust number, location and 
parameters of spots until photometric residua do not show any systematic 
or large-scale variability in phase. 

6. Relax a and q. Make a fit to the velocity curve. 

7. Fix a and q. Relax the remaining parameters, and make a fit to the 
photometric data. 

8. Go to (6). Repeat steps (6) - (8) until upon completing (7) the fit con- 
verges according to the same criterion as in (4). 

9. Repeat steps (5) - (8) for the third set of photometric data. 

10. Repeat steps (2) - (8) using another set of photometric data to obtain the 
base fit. 

11. Repeat steps (2) - (8) using the last set of photometric data to obtain the 
base fit. 
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An example result of fitting is shown in Fig. [21 with observational points omitted 
for clarity. In this case the base fit was obtained from the LCO data. The rms 
residual V magnitude of LCO, CTIO and SAAO data is equal to 0.007, 0.017 
and 0.016 mag, respectively; i.e., it is practically the same as the accuracy of 
photometric measurements reported in Section [2. 1.41 The rms residual B and I 
magnitudes amount to 0.014 and 0.017 mag, respectively. Note that all F-band 
light curves in Fig. [2] are asymmetric (the same holds for the B- and /-band 
light curves). The proximity effect is clearly visible in the "de-spotted" light 
curves (dotted lines in Fig. [2]), indicating a significant tidal distortion of the 
components of S892. 

The number of spots placed on model stars to match the photometric data 
varied from 1 to 5 per star in different fits. Both dark and bright spots were 
introduced, with |AToff |/Toff ranging from 0.01 to 0.06. In most cases the spot 



AV R 




0,05 - 



-0,05 - , . . . .■ . . . ■ :■ . 

0,05 - ^ ^ y" ' I ^ ^ ^ ^ ^ ' ^ ^ ' " ^ ' ^" I ^ ^ ' I ' ^ i ~ 

-0.6 -0.4 -0.2 0,2 0.4 0.6 

Phase (P = 0,451780 d) 

Figure 2: Upper panel: example fits to y-band light curves phased with P = 
0.451280 d (solid lines). Dotted lines: the same fits shown after the spots 
have been removed from both components of the system. From top to bottom: 
LCO, CTIO and SAAO. The zero-point on the vertical axis corresponds to 
y = 17.5 mag, and the curves are displaced from one another by 0.4 mag. Lower 
panels: residuals from fits with spots. From top to bottom: LCO, CTIO and 
SAAO. The base fit (see text for explanations) was obtained from LCO data. 
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Figure 3: Example modulations of y-band magnitude due to spots on pri- 
mary and secondary component of S892 (upper and lower panel, respectively). 
Dashed, dasli-dotted and solid lines: fits to LCO, CTIO and SAAO light curves. 

coverage was rather high - between 0.3 and 0.5 of the whole surface of the star. 
As the adopted configuration of spots is not unique we do not think it worthwhile 
to go into details, and we limit ourselves to just one example illustrated in Fig. [3] 
which shows periodic variations in the l^-band caused solely by the presence of 
spots (i.e., without the effect of ellipticity). The data for Fig. [3]are taken from 
models for which the base fit was obtained from the LCO photometry. 
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Phase (P = 0.451780 d) 

Figure 4: Example fit to UVES velocity data. The horizontal line shows the 
systemic velocity of —0.8 ±0.4 km s~^. 
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An example velocity fit for phases (f> < 0.9 is shown in Fig. 2] (bends on 
the synthetic curve are due to the Rossiter - McLaughlin effect). We decided 
to not include points with (j)> 0.95 despite their formal errors in Table [D being 
rather small because we suspected that heavy line-blending in the corresponding 
spectra would introduce unaccountable systematic deviations. Indeed, after a 
few trial fits with all velocity points taken into account we found that for 0>O.95 
it was impossible to get residual velocities smaller than 10 km s~^, whereas the 
rms residual velocity for (/)<0.9 was equal to only 2.1 km s~^. 



Table 3: Physical parameters of S892 



p 


[d] 


0.451780 


±0.000001 




a 


[Re] 


2.77 


±0.01 


/+0.013N 

l-o.oisJ 


e 


(fixed) 


0. 






i 


[degrees] 


85.36 


±0.28 


/+0.58\ 
\-0.24,) 


7 


[km s-i] 


-0.8 


±0.4 




mi 


[Mq] 


0.733 


±0.005 


/+0.011\ 
1-0.007.' 


7712 


[Me] 


0.662 


±0.005 


^+0.012\ 

l-o.oosJ 


Ri 


[Re] 


0.719 


±0.014 


^+0.030\ 
1-0.017J 


R2 


[Re] 


0.680 


±0.017 


^+0.031', 
1-0.014J 


Ti 


[K] (fixed) 


4200 






Ta 


[K] 


3940 


±20 


^+110x 

l-ioo.' 


Mbou 


[mag] 


6.89 


±0.04 


/-0.20-, 
V+0.15.' 


Afbol2 


[mag] 


7.29 


±0.04 


f-0.20\ 
(+0.16) 



3"^ column: rms deviations from mean values obtained in 9 fits to light and velocity curves. 
4**^ column above the horizontal line: largest deviations from the mean; below the horizontal 
line: variations corresponding to ±100 K change in Ti 

AU fits made according to the procedure described above produced 9 sets of 
system parameters whose averaged values are listed in Table [3] together with 
rms deviations from the nine-fit mean (since the " de-spotting" procedure may in 
principle introduce systematic errors, we also give the largest deviations from the 
mean). Obviously, both the mean values and the errors might change if we had 
more light curves. Unfortunately, a significant improvement of their accuracy 
would require at least ~10 additional light curves at epochs separated by at 
least several weeks. Note that temperatures - and, obviously, the bolometric 
magnitudes they yield - are not well constrained by the data (see Section |4]) . 
Their errors in column 3 of Table |3] are but formal results of the fits with 
Ti =4200 K; depending on the exact value of Ti they may vary as indicated in 
the last column. 

To check the influence of photometry errors on values given in Table [3] we 
removed spot effects from our best (LCO) and poorest (CTIO) data, and per- 
formed extensive (10000 points) Monte Carlo fitting with the help of the JK- 
TEBOP code (Southworth et al. 2007). For the LCO light curve the errors of 
i?i and i?2 were equal to 0.0021 Rq and 0.0026 Rq, respectively, while for the 
CTIO light curve they amounted to 0.0075 Rq and 0.0090 Rq. Since photom- 
etry errors Cph and "de-spotting" errors ads are independent, the total error 
may be approximated by the square root of the sum (o'pij±cr^s). Such an oper- 
ation causes the errors of Ri and i?2 to increase, respectively, by 14% and 13% 
compared to the values given in the 3'^'^ column of Table [31 

The radii of the components imply rotational velocities wi sini = 80.8 ± 
1.1 km s~^ and W2sini = 76.3 ±2.3 km s~^, which agree very well with ?;isini = 
79.2 ±5.9 km s^^ and t;2 sini = 76.2 ± 8.8 km s~^ derived from BF fits in Sec- 
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tion 12.21 thus indicating the consistency of our solution. Note also that the 
low systemic velocity of S892 adds credence to our assumption concerning its 
metallicity. 

4 Discussion 

Several authors, e.;;., Morales et al. (2008) and references therein, have found 
that radii of components of lower main sequence binary stars (LMSS) are sys- 
tematically and significantly larger than predicted by the theory. Our analysis 
of S892 confirms this finding. Incidentally, the radii of the components of this 
system agree very well with the empirical i?(M) relation 

R/Rq = 0.0324 + 0.9343M/Mq + 0.0374 {M/Mof (3) 

obtained by Bayless & Orosz (2006), which for Mi=O.733M0 and ==0.662Mq 
yields Ri — Q.THRq and R^ =O.667i?0. The updated M-R diagram for com- 
ponents of well-surveyed binary LMSS is shown in Fig. [5l 

As we already mentioned, based on our data it is impossible to find accurate 
values of Ti and T2. This is because we cannot derive absolute fluxes or effective 
temperatures from light and velocity data alone as we do not know the distance 
to S892 from independent measurements. In principle, we might try to match 
the spectra of this system to binary spectral templates following Becker et al. 
(2008). However, an additional problem arises here due to large and variable 
spot coverage of both components. As a consequence, the spectra of the compo- 
nents are intrinsically variable (see Fig. [S]), which puts any matching attempts 
in question. Some changes in Fig. [6] are due to orbital motion; however most 
of them (e.g., at 5460 < A [A] < 5480, 5530 < A [A] < 5590 and A[A]>5600) result 
from the intrinsic variability of the components. 

Another consequence of the vigorous activity of S892 is strong variability of 
the light curve, as illustrated in Fig. [T] One can see that at some orbital phases 
the luminosity of the system may vary by almost 10% within a year (note that 
apparent phase shifts in Fig. [7] are solely due to asymmetries of the light curves). 
Yet another aspect of the activity is strongly variable Ha emission from both 
components, with the equivalent width of the line changing within a few weeks 
from practically to ~ 7A per component (Fig. [8]) . H/3 is also seen in emission, 
with the equivalent width reaching up to ^3.5A. Neither Ha nor H/? are ever 
seen in absorption. We were not able to detect any secular or periodic intensity 
changes of either line, and we concluded that the hydrogen emission seems to 
be entirely chaotic. 

A system like this should be a fairly strong source of X-rays. Based on the 
relation between radius and ix/^boi found by Lopez-Morales (2007) we can 
expect the total X-ray luminosity of S892 to reach ~ 5 x 10^^ erg s^^. Unfortu- 
nately, because of the large distance of the system from the Sun the expected 
flux (~5 X 10~^^ erg cm~^ s""'^) is far too low to be registered in the ROSAT all- 
sky survey. We did not find relevant X-ray data in other catalogs or databases, 
either. An attempt to estimate the age of S892 based on the equivalent width 
of the Li I A6708 line was also futile - the line could not be detected. The qual- 
ity of our spectra is insufficient to draw any firm conclusion from this negative 
result - we can only say that the system almost certainly is not very young. 
Finally, we searched for proper motion data. S892 is listed in the UCAC3 Cat- 
alogue (Finch et al. 2009) as 3UC143-028414, with = -7.6 mas yr^^ and 
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Figure 5: Mass- Radius diagram for low-mass eclipsing binary stars with obser- 
vationally determined parameters (only stars with uncertainties below 3% both 
in mass and radius are included) . The components of S892 are indicated by dia- 
monds. The solid line is a theoretical 1.0 Gyr isochrone based on the solar com- 
position models of Baraffe et al. (1998). The data for NSVS 01031772 (marked 
with Ni^2), 2MASS J05162881+2607387 (marked with 2Mi^2), T-BooO-00080 
(marked with T-B) and T-Lyrl-01662 (marked with T-L) are taken, respec- 
tively, from Lopez-Morales et al. (2006), Bayless & Orosz (2006) and Fernandez 
et al. (2009). All remaining data (from left to right for CM Dra B, CM Dra A, 
CU Cue B, CU Cnc A, GU Boo B, YY Gem AB, GU Boo A, RXJ0239. 1-1028 B, 
RXJ0239. 1-1028 A and V818 Tau B) are taken from Ribas (2006). 

/i5 = — 3.8 mas yr~^, but large errors (6.4 mas yr^^ in both coordinates) make 
this information rather uninteresting. 

Based on composite CaH2, CaHa and TiOs molecular indices, Becker et al. 
(2008) derived = 3730 ± 100 K for the M = 0.66 M© primary component 
of 2MASS J01542930+0053266. While this resuh marginally agrees with our 
it is tempting to verify it using an additional evidence. Unfortunately, the 
indices used by Becker et al. (2008) are defined for A > 68OOA (Reid, Hawley 
& Gizis 1995), i.e., for wavelengths longer than the red limit of our spectra. 
Not being able to apply their method to our data, we compared their result to 
temperatures of similar systems reported in the literature. 

To our surprise, both stars in YY Gem and both stars in GU Boo, all being 
by -10% less massive than the primary of 2MASS J01542930-K0053266 were 
found to be by up to 350 K hotter (Torres & Ribas 2002; Lopez-Morales & Ribas 
2005). Active stars in this mass range are known to be cooler than nonactive 
ones (see Section [1]), but in all three systems the level of activity is similar (in 
principle, in 2MASS J01542930-h0053266 it should be even lower than in the 
remaining two because of a relatively long orbital period). Adopting T2 = 3730 K 
for our secondary and using the flux ratio obtained from our fits, we obtain for 
our 0.73 Mq primary Ti =3940 K - a value assigned by Lopez-Morales & Ribas 
(2005) to GU Boo A, whose mass is smaller by almost 20%. We conclude that 

■fNote that the larger error margin from the 4"" column of Table |3] should be applied here. 
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Figure 6: An illustration of the spectral variability of S892. From top to bottom: 
spectra obtained on HJD 2454746 ((/)=0.177), HJD 2454849 (0=0.695) and HJD 
2454744 (0 = 0.723). The flux scale is arbitrary; the top and bottom spectra are 
displaced vertically by 20 units. 
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Figure 7: Variability of the F-band light curve illustrated by differences between 
solid curves from Fig. H Dashed: CTIO-LCO; dash-dotted: SAAO-LCO; solid: 
CTIO-SAAO. 

3940 K from Table [3] is a fairly realistic value of T2 in the S892 system. Note 
that while all these discrepancies provide a good illustration of the notorious 
problem with temperature measurements, they luckily do not influence our main 
findings. 

Bolometric corrections BC\ = —1.00 mag and BC2 = —1.25 mag (Bessell, 
Castelli & Plez 1998) translate our AfboU and Mboi2 into Myi = 7.89 mag and 
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Figure 8: Ha emission in S892. The equivalent width of the line varies from 
to ^ 7A per component. The number in the upper right corner of each frame 
indicates the date of observation, in HJD-2454000. 



Mv2=8.54 mag, resulting in a combined absolute F-band magnitude of 7.41 mag 
at the phase of maximum light. NGC2204 is ^4 kpc distant (Kassis et al. 
1997) and only weakly reddened, with E{B-V)r^0.l (Schlegel et al. 1998) and 
Ay ~ 0.32. Thus, the revised distance to S892 is about 1050 pc if the whole 
absorbing material is located behind the system, and about 900 pc if it is located 
in front of it. 



5 Conclusions 

We discovered and analyzed a detached eclipsing binary composed of two late 
K dwarfs. Based on light and velocity curves we derived orbital period of 
0.451780±0.000001 d, semimajor axis a=2.77±0.01 Rq and inclination i=85?36 
±0?28. The masses of the components (Mi=0.733±0.005 Mq- M2=0.662±0.005 
Mq) and their radii (iii = 0.72 ± 0.01 Rq\ i?2 =0.68 ±0.02 Rq) are consistent 
with the empirical mass-radius relationship established recently for LMSS in 
binary systems by Bayless & Orosz (2006): both stars are larger and cooler 
compared to theoretical models. 

The discrepancy between theory and observations was originally thought to 
arise from inaccuracies in the equation of state and/or opacity. It turned out, 
however, that models based on the same physics were entirely compatible with 
observed properties of single stars or stars in wide binaries (Lopez-Morales et al. 
2006). Another approach to the problem was suggested by a positive correlation 
between radius and activity level found for a few members of close binaries by 
Lopez-Morales (2007). It was realized that the basic difference between single 
stars and components of short-period binary systems was high rotational veloc- 
ities of the latter, promoting high levels of magnetic activity. Strong magnetic 
fields can inflate the star due to large spot coverage (thus lowering the effective 
temperature) and/or partial inhibition of convection (which has the same effect 
as increased opacity; see Chabrier et al. for details). 

Both components of our binary are very active, showing high spot-coverage, 
unstable light curves, and strong variability of the emission in hydrogen lines. 
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Thus, our results support the increasing evidence that the observed inflation of 
radii of K and M dwarfs is related to high levels of magnetic activity. 
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